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We examine the use of string diagrams and the mathematics of category theory in 
the description of quantum states by tensor networks. This approach lead to a unifi- 
cation of several ideas, as well as several results and methods that have not previously 
appeared in either side of the literature. Our approach enabled the development of a 
tensor network framework allowing a solution to the quantum decomposition prob- 
lem which has several appealing features. Specifically, given an n-body quantum 
state we present a new and general method to factor l^) into a tensor network 
of clearly defined building blocks. We use the solution to expose a previously un- 
known and large class of quantum states which we prove can be sampled efficiently 
and exactly. This general framework of categorical tensor network states, where a 
combination of generic and algebraically defined tensors appear, enhances the theory 
of tensor network states. 
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I. INTRODUCTION 

Tensor network states have recently emerged from Quantum Information Science as a gen- 
eral method to simulate quantum systems using classical computers. By utilizing quantum 
information concepts such as entanglement and condensed matter concepts like renormaliza- 
tion, several novel algorithms, based on tensor network states (TNS), have been developed 
which have overcome many pre-existing limitations. These and other related methods have 
been used to perform highly accurate calculations on a broad class of strongly-correlated 
systems and have attracted significant interest from several research communities concerned 
with computer simulations of physical systems. 

In this work we develop a tool set and corresponding framework to enhance the range 
of methods currently used to address problems in many-body physics. In this categorical 
network model of quantum states, each of the internal components that form the build- 
ing blocks of the network are defined in terms of their mathematical properties, and these 
properties are given in terms of equations which have a graphical interpretation. In this 
way, diagrammatic methods from modern algebra and category theory |l| can be combined 
with graphical methods currently used in tensor network descriptions of many-body physics. 
Moreover our results indicate that it may be be advantageous to use "categorical compo- 
nents" within tensor networks, whose algebraic properties can permit a broader means of 
rewiring networks, and potentially reveal new types of contractible tensor networks. Our 
results include defining a new graphical calculus on tensors, and exposing their key proper- 
ties in a "tensor tool box" . We use these tensors to present one solution to the problem of 
factoring any given quantum state into a tensor network. The graphical properties of this 
solution then enabled us to expose a wide class of tensor network states that can be sampled 
in the computational basis efficiently and exactly. Here we list several novel contributions 
of the present paper. 

(i) We introduce a universal class of tensors to the quantum theory and tensor networks 
literature. We cast several properties and simplification rules applicable for classical 
networks of this type, into the language of tensor networks. (This fixed collection of 
tensors and the corresponding Boolean algebra framework we introduce provides the 
potential for a new tool to construct networks of relevance for problems in quantum 
information science and condensed matter physics. Recent work has used the collection 
of tensors we introduce here (arXiv version) together with their rewrite rules to solve 
models in lattice gauge theory [2j and to study correlator product states (3)]). 

(ii) We define a new class of quantum states, quantum Boolean states, expressed as: = 
Y2f( x )\ x ) where / is a switching function and the sum is over all n-long bit strings. 
A quantum state is called Boolean iff it can be written in a local basis with amplitude 
coefficients taking only binary values or 1. Examples of states in this class include 
states such as W and GHZ -states. We develop tools which aid in the study of this 
class of states. 

(iii) We prove that a tensor network representing a Boolean quantum state is determined 
from the classical network description of the corresponding function. The quantum 
tensor networks are found by the following method: we let each classical gate act on a 
linear space and replace the composition of functions, with the contraction of tensors. 
This technique is detailed in the present work and has very recently been used in [2}. 

(iv) We present a proof which shows a new universal constructive decomposition of any 
quantum state, into a network built in terms of tensors in the tool box we introduce. 
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(v) We present a new and large class of quantum states which we prove to be efficient to 
sample. A subclass of this class of states includes the widely studied class of correlator 
product states. This connection (arXiv version) was recently studied in jjjj]. Exploring 
this class in further detail is left to future studies. 

(vi) Although Boolean algebra has long been used in both classical and quantum computer 
science, we seem to be the first to tailor its use to the problem of describing quantum 
states by tensor networks. This has resulted in new proof techniques and also enables 
one to use methods from the well developed graphical language of classical circuits 
inside the domain of tensor networks. 

To explain the main motivation which prompted us to study tensor network states, let us 
recall the success of established numerical simulation methods, such as the density matrix 
renormalization group (DMRG) [U-0] which is based on an elegant class of tensor networks 
called Matrix Product States (MPS) Q- For more than 15 years DMRG has been a key 
method for studying the stationary properties of strongly-correlated ID quantum systems in 
regimes far beyond those which can be described with perturbative or mean-field techniques. 
Exploiting the tensor network structure of MPS has lead to explicit algorithms, such as 



the Time- Evolving Block Decimation (TEBD) method |8Hl3|. for computing the real-time 
dynamics of ID quantum systems. Accurate calculations of out-of-equilibrium properties 



prope 

has proven extremely useful for describing various condensed matter systems [14KL6 . as 
well as transport phenomena in ultra-cold atoms in optical lattices [l7- 2 . Additionally the 
TEBD method has recently been successfully adapted to the simulation of stochastic classical 
systems [20], as well as for simulating operators in the Heisenberg picture 0. Despite these 
successes, limitations remain in the size, dimensionality, and classes of Hamiltonians that 
can be simulated with MPS based methods. To overcome these restrictions, several new 
algorithms have been proposed which are based on different types of tensor network states. 



Specifically: Projected Entangled Pair States (PEPS) [2ll423l| which directly generalize the 
MPS structure to higher dimensions, and the Multi-scale Entanglement Renormalization 
Ansatz (MERA) [24142 7|) which instead utilizes an intuitive hierarchical structure. 



Category theory is often used as a unifying language for mathematics []J and in more 



recent times has been used to formulate physical theories [28( . One of the strong points of 
categorical modeling is that it comes e quipped with many types of intuitive graphical calculi. 
(We mention the coherence results [j], 29, 30] and as a matter of convenience, make use of 



f-compactness. The graphical calculus of categories formally extends to a rigorous tool. 
See for instance, Selinger's survey of graphical languages for monoidal categories outlining 
the categories describing quantum theory (30[.) The graphical theory behind the types of 
diagrams we consider here dates to the work of Lafont |3lL 32| who built directly on the 
earlier work by Penrose [33| . 



A motivating reason for connecting category theory to tensor networks is that, increas- 
ingly, both existing and newly developed tensor network algorithms are most easily expressed 
in terms of informal graphical depictions. This graphical approach can now be comple- 
mented and enhanced by exploiting the long existing language of category theory [l, 28|, 34 



This immediately enables the application of many established techniques allowing for both 
a "zoomed out" description exposing known high-level structures, but also enables new 
"zoomed in" descriptions, exposing "hidden" algebraic structures that have not previously 
been considered. 

We will illustrate our categorical approach to tensor networks by focusing on tensor net- 
works constructed from familiar components, namely Boolean logic gates (and multi-valued 
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logic gates in the case of qudits), applied to this unfamiliar context. To accomplish this 
goal, we build on ideas across several fields. This includes extending the work by Lafont [35[ 
which was aimed at providing an algebraic theory for classical circuits. (Lafont 's work is 



related to the more recent work on proof theory by Guiraud [36j, and is a different direction 
from other work on applying category theory to classical networks appearing in (37j].) The 
use of symmetric monoidal categories tightens this approach and removes some redundancy 
in Lafont 's graphical lemmas. The application of these results to tensor networks introduces 
several novel features. The first feature is that once conventional logic circuits are formulated 
as tensor networks they can be distorted into atemporal configurations since the indices (or 
legs) of tensors can be bent around arbitrarily. This permits a very compact tensor network 
representation of a large interesting class of Boolean states such as GHZ-states, W-states and 
symmetric states [38|], using exclusively Boolean gate tensors. A second feature is that once 
expressed as tensors the corresponding classical logic circuits act on complex valued inputs 
and outputs, as opposed to just binary values. By permitting arbitrary single-qubit states 
(general rank-1 tensors) at the output of tensor networks, which are otherwise composed 
of only switching functions, we arrive at a broader class of generalized Boolean states. We 
prove in Theorem [23] that this class of states provides an explicit construction method for 
factoring any given quantum state into a tensor network. As expected, the cost of this 
exhaust iveness is that the resulting network is, in general, neither efficient in description 
or in contraction. However, by limiting both the gate count and number of the switching 
functions comprising the tensor networks to be polynomial in the system size, we obtain a 
class of states, which we call Generalized Polynomial Boolean States (GPBS see Definition 
|2"51) . that can be sampled in the computational basis efficiently and exactly (Theorem I2"9j) . 

a. Manuscript Structure. Next in Section [III we quickly review the key concepts in- 
troduced in this paper before going into detail in the remaining sections. We continue in 
Section II III by defining the network building blocks: this includes defining some new rank- 
3 tensors such as the quantum AND-state in Equation (jUJ). We then consider how these 
components interact in Section IIVI This is done in terms of algebraic structures, such as 
Bialgebras (Section HVBj) and Hopf-algebras (Section IIVB II) which are well known to have 
a purely diagrammatic interpretation. With these definitions in place, in Section [V] we apply 
this framework to tensor network theory. As an illustrative example, we consider the W- 
and GHZ-states using our formalism. Specifically, we consider a particular categorical tensor 
network for many-body W-states in Section [V] A proof of our decomposition theorem for 
quantum states is given in Section IVH In conclusion, we mention some future directions for 
work in Section IVHI We have included Appendix [A] on algebras defined on quantum states 
and Appendix [B] on the Boolean XOR-algebra. 

b. Background Reading. The results appearing in this work were found by tailoring 
several powerful techniques from modern mathematics: category theory, algebra and co- 
algebra and applicable results from classical network theory and graphical calculus. Tensor 
network states are covered in the reviews For general background on category 



theory see [l[ . For background on Boolean algebra, discrete set functions and circuit theory 
see 41] and see 42f for background on pseudo Boolean functions and for multi- valued logic 
see 43]. For background on quantum circuits and quantum computing concepts see jH, 45| 
and for background on the theory of entanglement see For the current capabilities of 



the existing graphical language of tensor network states see e.g. [471 . |48j and for work on 



using ideas related to tensor networks for state preparation of physical systems see [48H50 . 
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II. RESULTS OVERVIEW 



In the present Section, we informally review our main results. The idea of translating 
any given quantum state or operator into a representation in terms of a connected network 
of algebraically defined components is explained next in Section III Al with the correspond- 
ing algebraic definitions of these network components over viewed in Section III Bl Boolean 
quantum logic tensors are then introduced in Section III CI We summarize our main results 
in Section Hi Dl 



A. Tensor network representations of quantum states 

A qudit is a d- level generalization of a qubit. In physics a quantum state of n-qudits has an 
exact representation as a rank-n tensor with each of the open legs corresponding to a physical 
degree of freedom, such as a spin with {d— l)/2 energy levels. Such a representation, shown 
in Figure d^a) is manifestly inefficient since it will have a number of complex components 
which grows exponentially with n. The purpose of tensor network states is to decompose 
this type of structureless rank-n tensor into a network of tensors whose rank is bounded. 
There are now a number of ways to describe strongly-correlated quantum lattice systems 
as tensor-networks. As mentioned in the introduction, these include MPS 0, SI, Hi 



PEPS |39|,|53| and MERA j25|,[54. For MPS and PEPS, shown in Figures QJb) and (c), the 



resulting network of tensors follows the geometry of the underlying physical system, e.g., a 



ID chain and 2D grid, respectively. Alternatively a Tensor Tree Network (TTN) [55|, [56 
can be employed which has a hierarchical structure where only the bottom layer has open 
physical legs, as shown in Figure [T^d) for a ID system and Figure (He) for a 2D one. (Each 
tensor in these networks is otherwise unconstrained, although enforcing some constraints, 
such as orthogonality, has numerical advantages.) For MERA the network is similar to 
a TTN, as seen in Figure [T]^f) for ID, but is instead comprised of alternating layers of 
rank-4 unitary and rank-3 isometric tensors. The central problem faced by all types of 
tensor networks is that the resulting tensor network for the quantity (ip\(0\ip)) , where O 
is some product operator, needs to be efficiently contractible if any physical results, e.g., 
expectation values, correlations or probabilities, are to be computed. For MPS and TTN 
efficient exact contractibility follows from the ID chain or tree-like geometry, while for MERA 
it follows from its peculiar causal cone structure resulting from the constraints imposed on 
the tensors [25| . For PEPS, however, exact contraction is not efficient in general, but can 



be rendered efficient if approximations are made [39|, [53 . 

In our approach we define a categorical tensor network state (CTNS) generally as any 
TNS which contains some algebraically constrained tensors along with possible generic ones. 
Indeed, when recast, certain widely used classes of TNS can be readily exposed as examples 
of CTNS. Specifically, variants of PEPS have been proposed called string-bond states [571 ]. 
Although these string-bond states, like PEPS in general, are not efficiently contractible, they 
are efficient to sample. By this we mean that for these special cases of PEPS any given am- 
plitude of the resulting state (for a fixed computational basis state) can be extracted exactly 
and efficiently, in contrast to generic PEPS. This permits variational quantum Monte-Carlo 
calculations to be performed on string-bond states where the energy of the state is stochas- 
tically minimized [57] . This remarkable property follows directly from the use of a tensor, 
called the COPY-dot, which will form one of several tensors in the fixed toolbox considered 
in great detail later. As its name suggests, the COPY-dot duplicates inputs states in the 



6 



(a) 



CD 




(e) 



(b) -^pp^pppp^- (d) 




(f) 




FIG. 1. (a) A generic quantum state \tp) for n degrees of freedom represented as a tensor with n 
open legs, (b) A comb-like MPS tensor network for a ID chain system 5l|, |52|. (b) A grid- like 
PEPS tensor network for a 2D lattice system 39,153]. (d) A TTN for a ID chain system where only 
the bottom layer of tensors possess open physical legs 55j, [56|]. (e) A TTN for a 2D lattice system, 
(f) A hierarchically structured MERA network for a ID chain system possessing unitaries (rank-4 
tensors) and isometries (rank-3 tensors) 25j, |5J|. This tensor network can also be generalized to a 
2D lattice (not shown). 



computational basis, and thus with these inputs breaks up into disconnected components, as 
depicted in Figure [2]^a). By using the COPY-dot as the "glue" for connecting up a TNS, the 
ability to sample the state efficiently is guaranteed so long as the individual parts connected 
are themselves contractible. The generality and applicability of this trick can be seen by 
examining the structure of string-bond states, as well as other types of similar states like 
entangled-plaquette-states [58f and correlator-product states j59[, shown in Figure Ws)-{e). 
A long-term aim of this work is that by presenting our toolbox of tensors, entirely new 
classes of CTNS with similarly desirable contractibility properties can be devised. Indeed 
we have a useful result in this direction (Theorem 129]) by introducing a new class of states, 
which can also be sampled exactly and efficiently. 



On an interesting historical note, to the best of our knowledge, a graphical interpretation 
of tensors was first pointed out in 33[ . A graphical language for describing the manipulations 
and steps of tensor network based algorithms has become widely used. By introducing new 
tensors, and by considering their graphical properties, an aim of our work is to extend the 
existing methods of the diagrammatic methods used. The graphical properties of tensors 
are defined succinctly via so-called string diagrams. As an exemplary illustration, we will 
consider in great detail CTNS which are composed entirely from a tensor toolbox built from 
classical Boolean logic gates. By invoking known theorems asserting the universality of 
multi- valued logic [43| (also called d-state switching), our methods can be readily applied to 
tensors of any finite dimension. Our approach provides not only an example of the use of 
well known gates in an unfamiliar context but also illustrates the potential power of having 
"algebraic components" within a tensor network. The next two sections highlight some of 
the properties of this toolbox of tensors and reviews our main result showing a new quantum 
state decomposition using a subset of them. The full details of this work then follow from 
Section 0TO onwards. 
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FIG. 2. (a) One of the simplest tensors, called the diagonal in category theory, the COPY-gate 
or the COPY-dot in circuits, copies computational basis states \x) where x = 0, 1 for qubits and 
x = 0, 1, ...,d — 1 for qudits. The tensor subsequently breaks up into disconnected states, (b) A 
generic PEPS in which we expose a single generic rank-5 tensor. This tensor network can neither 
be contracted nor sampled exactly and efficiently. However, if the tensor has internal structure 
exploiting the COPY-dot then efficient sampling becomes possible, (c) The tensor breaks up into 
a vertical and a horizontal rank-3 tensor joined by the COPY-dot. Upon sampling computational 
basis states the resulting contraction reduces to many isolated MPS, each of which are exactly 
contractible, for each row and column of the lattice. This type of state is known as a string-bond 
state and can be readily generalized [57|. (d) An even simpler case is to break the tensor up into 
four rank-2 tensors joined by a COPY-dot forming a co-called correlator-product state [59[]. (e) 
Finally, outside the PEPS class, there are entangled plaquette states [3] which join up overlapping 
tensors (in this case rank-4 ones describing a 2 x 2 plaquette) for each plaquette. Efficient sampling 
is again possible due to the COPY-dot. 



B. Network components fully defined by diagrammatic laws 

We will now review the set of tensors that form our universal building blocks. To get an 
idea of how the tensor calculus will work, consider Figure [31 which forms a presentation of 
the linear fragment of the Boolean calculus (35[): that is, the calculus of Boolean algebra we 
represent on quantum states, restricted to the building blocks that can be used to generate 
linear Boolean functions. 

To recover the full Boolean-calculus, we must consider a non-linear Boolean gate: we use 



the AND-gate. Figure |3] together with Figure H] form a full presentation of the calculus [35 
The origin and consequences of these relations will be considered in full detail in Section [TTTl 
The presentations in Figure [3] together with Figure H] represent a complete set of defining 
equations [6(| • The results we report and the introduction of this new picture calculus into 
physics has already attracted significant interest and provided a new research direction in 
categorical quantum mechanics [6ll . 62 



Proceeding axiomatically we need to add a bit more to the presentation of the Boolean 
calculus to represent operators and quantum states. This is because e.g. all the diagrams 
in Figure [3] and H] are read from the top of the page to the bottom. Our network model 
of quantum states requires that we are able to turn maps upside down, e.g. transposition. 
This additional flexibility comes from an added ability to bend wires. We can hence define 



s 




FIG. 3. Read top to bottom. A presentation of the linear fragment of the Boolean calculus. The 
plus (©) dots are XOR and the black (•) dots represent COPY. The details of (a)-(g) will be given 
in Sections IIIII and IIVI For instance, (d) represents the bialgebra law and (g) the Hopf-law (in the 
case of qubits x © x = 0, in higher dimensions the units (+| becomes (0| + (1| + • • • + {d — 1|). 




FIG. 4. Diagrams read top to bottom. A presentation of the Boolean-calculus with Figure [3l The 
details of (a)-(g) will be given in Sections IIIII and HVl For instance, (h) represents distributivity of 
AND(A) over XOR (©), and (d) shows that x Ax = x. 

transposition graphically (see Figure [181 (d)). 

The way forward is to add what category theory refers to as compact structures (see 
Section HV CI for details). These compact structures are given diagrammatically as 
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CD 




1 f 

and as will be explored in Section TlV CI these two structures allow us to formally bend wires 
and to define the transpose of a linear map/state, and provide a formal way to reshape a 
matrix. We understand (a) above as a cup, given as the generalized Bell-state Si=o 1^) 
and (b) above as the so-called cap, Bell-costate Si=o or e Jf ec ^- (Normalization factors 
omitted: without loss of generality, we will often omit global scale factors (tensor networks 
with no open legs). This is done for ease of presentation. We note that for Hilbert space H 
there is a natural isomorphism C ® H = H = 7i <S> C.) 

Compact structures provide a formal way to bend wires — indeed, we can now connect a 
diagram represented with an operator with spectral decomposition SiANK^I bend all the 
open wires (or legs) towards the same direction and it then can be thought of as representing 
a state (Si&K)K) where overbar is complex conjugation), bend them the other way and 
it then can be thought of as representing a measurement outcome (Si /3i(i\(i\), that is 
an effect. (The isomorphism SiA(*|(*| — SiAK)(*l — Si f° r a rea l valued basis 

becomes Si A^'K^I — Si — Si which amounts to nipping a bra to a ket and 
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00) + |01) + |10) 



FIG. 5. Example of the Boolean quantum AND-state or tensor. In (a) the network is run backwards 
(post-selected) to (1| resulting in the product state |11). In (b) the tensor is post-selected to (0| 
resulting in the entangled state |00) + |01) + 1 10) . 

vise versa.) One can also connect inputs to outputs, contracting indices and creating larger 
and larger networks. With these ingredients in place let us now consider the new class of 
Boolean quantum states. 

C. Boolean and multi- valued tensor network states 

To illustrate the idea of defining Boolean and multi-valued logic gates as tensors, consider 
Figure [5] which depicts a simple but key network building block: the use of the so-called 
"quantum logic AND-tensor" which we define in Section [HI D[ This is a representation of 
the familiar Boolean operation in the bit pattern of a tri-qubit quantum state as 

|^and)= ki) «) |ar 2 > ® A ac 2 ) = |000> + |010> + |100> + |111) 

x 1 ,x 2 &{0,l} 

and hence the truth table of a function is encoded in the bit pattern of the superposition 
state. This utilizes a linear representation of Boolean gates on quantum states as opposed 
to the typical direct sum representation common in Boolean algebra. 

In this work we are particularly concerned with network constructions as a means to 
study many-body quantum states by tensor networks. First, we can compose AND-states 
(by connecting wires and hence contracting tensor indices) — together with NOT-gates, this 
enables one to create the class of Boolean states in Equation (j2j). That is, one will realize 
a network that outputs logical-one (represented here as |1)) any time the input qubits 
represent a desired term in a quantum state (e.g. create a function that outputs logical-one 
on designated inputs 1 00) , 1 01) and 1 10) and zero otherwise as shown in Figure \5§. We then 
insert a |1) at the network output. This procedure recovers the desired Boolean state as 
illustrated in Figure E^a) with the resulting state appearing in Equation ([T]). 

(1) 

The network representing the circuit is read backwards from output to input. Alternatively 
the full class of Boolean states is defined as: 

Definition 1 (Boolean many-body qudit states). We define the class of Boolean states as 
those states which can be expressed up to a global scalar factor in the form (j2J) 



) y \xi,x 2 , ■■.,x n )\f(x 1 ,x 2 , -;X n )) 

xi,X2,—,x n e{0,l,...,d— 1} 



(2) 
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where / : — > Z^ is a <i-switching function and the sum is taken to be over all variables Xj 
taking and 1 for qubits and 0, 1, ...d — 1 in the case of d-level qudits (see Figure |6] (a)). 




FIG. 6. A general multi-valued qudit state based on a d-switching function / can either be formed 
as (a) by inputting a logical-one at the output of the circuit as described by Equation (prj or (b) 
by bending the output of the circuit around to form an input as in Equation ([2]) . For multi- valued 
qudit states, the boolean function / : Z?> — >■ Z2 becomes a multi-valued qudit function / : — > Z^. 
The network (a) is then post-selected to oo|0) +ai|l) + • • -ad-i\d— 1) where Vz,Oj = 0/1. 

Examples of Boolean states include the familiar GHZ-state 1 00 • • • 0) + 1 11 • • • 1) which on 
qudits in dimension d becomes 

d-l 

|GHZ d ) = £ = |0)|0)|0) + |1)|1)|1) + • • • + |d- l)\d- l)\d - 1) (3) 

i=0 

as well as the W-state |00 ■ • • 1) + 101 ■ ■ ■ 0) + ■ ■ ■ + 1 10 ■ ■ ■ 0) which again on qudits becomes (in 
Equation (j3J) the operator X\m) = \m+ l(mod d)) is one way to define negation in higher 
dimensions. The subscript labels the ket (labeled 1,2 or 3 from left to right) the operator 
acts on i times.) 

d-l 3 

l W <*> : = EE W|0)|0)|0) = |0)|0)|1) + |0)|1)|0) + |1)|0)|0) + |0)|0)|2) + |0)|2)|0) + |2)|0)|0) + • • • 
t=i j=i 

(4) 

••• + |0)|0)|d-l) + |0)|d-l)|0) + \d- 1)|0)|0) 

What is clear from this definition is that Boolean states are always composed of equal 
superpositions of sets of computational basis states, as the allowed scalars take binary values, 
0,1. Our main result is that, despite this apparent limitation, tensor networks composed 
only of Boolean components can nonetheless describe any quantum state. To do this we 
require a minor extension to include superposition input/output states, e.g. rank-1 tensors 
of the form 1 0) + /3i 1 1) + • • • + fid- i I d — 1 ) . This gives a universal class of generalized Boolean 
tensor networks which subsumes the important subclass of Boolean states. This class is then 
shown to form a nascent example of the exhaustiveness of CTNS and to give rise to a wide 
class of quantum states that we show are exactly and efficiently sampled. 

D. Putting it all together: connecting the dots 

The key point to this result is that the introduction of Boolean logic gate tensors into 
the tensor network context allows the seminal logic gate universality results from classical 
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network theory to be applied in the setting of tensor network states. By extending this result 
we can construct a solution to the related quantum problem — that is, the decomposition 
or factorization of any quantum state into a CTNS. Thus our main result is captured by 
the following statement (see Theorem 1251) . 

Result (Translating quantum states into categorical tensor networks). Given quantum 
state Theorem [23] asserts a constructive method to factor into a CTNS constructed 
from rank-3, rank-2 tensors taken solely from the fixed set in the presentation from Figure [3] 
and H] together with arbitrary rank-1 tensors. 

This example then demonstrates the exhaust iveness of the most extreme case of the 
CTNS approach, where almost all tensors are chosen from a small fixed set of tensors with 
precisely defined algebraic properties. Importantly, in Theorem [29] the form of this general 
construction is limited in such a way as to provide a new class of states which can be exactly 
and efficiently sampled. 



III. CONSTITUENT NETWORK COMPONENTS: A TENSOR TOOL BOX 

Any vector space V has a dual V*: this is the space of linear functions / from V to the 
ground field C, that is / : V — >■ C. This defines the dual uniquely. We must however fix a 
basis to identify the vector space V with its dual. Given a basis, any basis vector \i) in V 
gives rise to a basis vector (j\ in V* defined by = b\ (Kronecker's delta). This defines 
an isomorphism V — > V* sending |i) to (i\ and allowing us to identify V with V*. In what 
follows, we will fix a particular arbitrarily chosen basis (called the computational basis in 
quantum information science). We will now concentrate on Boolean building blocks that 
are used in our construction. 



A. COPY-tensors: the "diagonal' 



The copy operation arises in digital circuits |41|, [63j and more generally, in the context 
of category theory and Algebra, where it is called a diagonal in cartesian categories. The 
operation is readily defined in any finite dimension as 



d-l 

def 



A d ^]>»<i| (5) 



i=0 



As |0) and |1) are eigenstates of a z , we might give A the alternative name of Z-copy. In 
the case of qubits COPY is succinctly presented by considering the map A that copies 
cr^-eigenstates: 

>}h+|00> 



a : c -> e 



"2 .. 



|i>-Hii> 



This map can be written in operator form as A : |00)(0| + |11)(1| and under cup/cap induced 
duality (on the right bra) this state becomes a GHZ-state as IV'ghz) = 1 000) + 1 1 1 1 ) = 
1 00} (0 1 + |11)(1|. The standard properties of COPY are given diagrammatically in Figure [7] 
and a list of its relevant mathematical properties are found in Table [B 



12 



(a 




FIG. 7. Salient diagrammatic properties of the COPY-dot. (a) Full-symmetry, (b) Copy points, 
e.g. \x) i— > \xx) for x = 0, 1 for qubits and x = 0, 1, d — 1 for d dimensional qudits. (c) The unit 
- in this case the unit corresponds to deletion, or a map to the terminal object which is given 
as (+| = f (0[ + (1| for qubits and (+| = f (0| + (1| + • • • + (d — 1| for d dimensional qudits (the bi- 
direction of time is explained later by considering co-diagonals in Section Till E|) . (d) Co-interaction 
with the unit creates a Bell state Yli=o K*)- This and the corresponding dual under the dagger 
form the compact structures of the f-category of quantum theory. 

Remark 2 (The COPY-gate from CNOT). The CNOJ-gate is defined as |0)(0| 1 <g) 1 2 + 
® erf. We will set the input that the target acts on to |0) then calculate CNOT(l! © 
|0) 2 ) = |0)(0| 1 © |0) 2 + |1)(1| 1 © |1) 2 . We have hence defined the desired COPY map copying 
states from the Hilbert space with label 1 (subscript) to the joint Hilbert space labeled 1 and 
2. 



B. XOR-tensors: the "addition" 



The XOR-gate implements exclusive disjunction or addition (mod 2 for qubits) and is 

We note that for multi-valued logic a modulo subtraction 
. By what could be called "dot-duality" the XOR-gate is 



denoted by the symbol © |64|, [65 
gate can also be denned as in 



66 



simply a Hadamard transform of the COPY-gate, appropriately applied to all of the dots legs. 



z i27rab / d \a)(b\ 



H * 



(We denote the discrete Fourier transform gate by Hu := ^Xla6e{oi <z-i} 
where d = dim'H is the dimension of the Hilbert space the gate acts in. We can see that 
H T = H, and that in a qubit system H coincides with the one-qubit Hadamard gate 66].) 
This can be captured diagrammatically in the slightly different form: 





To define the gate on the computational basis, we consider f{x\, x 2 ) = X\ © x 2 then / = 
corresponds to (xi,x 2 ) G {(0, 0), (1, 1)} and / = 1 corresponds to (xi,x 2 ) G {(1, 0), (0, 1)}, 
where the truth table for XOR follows 



Xi 


x 2 


f(x 1 ,x 2 ) = Xi @x 2 














1 


1 


1 





1 


1 


1 






Under cap/cap induced duality, the state defined by XOR is given as 



def 



E 



x 1 )\x 2 )\f{x 1 ,x 2 )) = |000) + |110) + |011) + |101) 



(6) 



13 



which is in the GHZ-class by LOCC equivalence viz. \if> @ ) = H g> H © H(|000) + |111)). The 
operation of XOR is summarized in Table UTI Since the XOR-gate is related to the COPY-gate 
by a change of basis, its diagrammatic laws have the same structure as those illustrated in 
Figure [3 The gate acting backwards (co-XOR) is defined on a basis as follows: 

^C®C::< l °\^ m + \ n \ or equivalently 
^ |1) ^ |10) + |01) 




C. Generating the affine class of networks 

Thus far we have presented the XOR- and COPY- gates. This system allows us to create 
the linear class of Boolean functions. As explained in the present subsection, this class can 
be extended to to the affine class by introducing either a gate that acts like an inverter, or by 
appending a constant |1) into our system. This constant will allow us to use the XOR-gate 
to create an inverter. 

A complemented Boolean variable is a Boolean variable that appears in negated form, 
that is -ix or written equivalently as x. Negation of a Boolean variable x can be expressed 
as the XOR of the variable with constant 1 as x = 1 © x. Whereas uncomplimneted Boolean 
variables are Boolean variables that do not appear in negated form (e.g. negation is not 
allowed). Linear Boolean functions contain terms with uncomplemented Boolean variables 
that appear individually (e.g. variable products are not allowed such as X\X 2 and higher 
orders etc., see Section |B]). Linear Boolean functions take the general form 

f(x 1 ,x 2 , x n ) = CiXi © c 2 x 2 © ... © c n x n (7) 

where the vector (ci, c 2 , c n ) uniquely determines the function. The affine Boolean func- 
tions take the same general form as linear functions. However, functions in the affine class 
allows variables to appear in both complemented and uncomplemented form. Affine Boolean 
functions take the general form 

f(xi, X 2 , X n ) = C © CiXi © C 2 X 2 © ... © C n X n (8) 

where c = 1 gives functions outside the linear class. From the identities, 1 © 1 = and 
© x = x we require the introduction of only one constant (co), see Appendix IB1 

Together, XOR and COPY are not universal for classical circuits. When used together, 
XOR- and COPY-gates compose to create networks representing the class of linear circuits. 
The affine circuits are generated by considering the constant |1). The state |1) is indeed 
copied by the black dot. However, our axiomatization (Figure [3]) proceeds through consid- 
ering the XOR- and COPY-gates together with |+), the unit for COPY and |0) the unit for 
XOR. It is by appending the constant |1) into the formal system (Figure [3]) that the affine 
class of circuits can be realized. 

Remark 3 (Affine functions correspond to a basis). Each affine function is labeled by a 
corresponding bit pattern. This can be thought of as labeling the computational basis, as 
states of the form |{0, 1}™) are in correspondence with polynomials in algebraic normal form 
(see Appendix lB]) . 
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D. Quantum AND-state tensors: Boolean universality 

The proceeding sections have introduced enough machinery to generate the linear and 
affine classes of classical circuits. These classes are not universal. To recover a universal 
system we must add a non-linear Boolean gate. We do this by representing the AND gate 
as a tensor. The unit for this gate is (1| and so can be used to elevate the linear fragment 
to the affine class. 




The AND gate (that is, A) implements logical conjunction [41|, |63|. Using again "dot- 
duality", the AND-gate relates to the OR-gate via De Morgan's law. This can be captured 
diagrammatically as 

To define the gate on the computational basis, we consider /(xi,x 2 ) = X\ A x 2 which we 
write in short hand as X\X 2 . Here / = corresponds to (xi, x 2 ) G {(0,0), (0, 1), (1,0)} and 
/ = 1 corresponds to (xi,x 2 ) = (1, 1). 

Under cap/cap induced duality, the state defined by AND is given as 

|Va) = Yl \xi)\x 2 )\f(x llX2 )) = |000) + |010) + |010) + |111) (9) 

xi,x 2 e{o,i} 

The key diagrammatic properties of AND are presented in Figure Eland the gate is summa- 
rized in Table IITT1 The gate acting backwards (co-AND) is defined on a basis as follows: 

A:C^CWJ!°^! 00 ! + |01> + |10> or (! + ^' + + \ , v , y 

|i> i— s- |n> 1 |-}i-> |oo> + |oi> + |io> -|ii) 




FIG. 8. Salient diagrammatic properties of the AND-tensor. (a) Input-symmetry, (b) Existence of 
a zero or fixed-point, (c) The unit |1). (d) Co-interaction with the unit creates a product-state. 
Note that the gate forms a valid quantum operation when run backwards as in (d). 



Example 4 (AND-states from Toffoli-gates) . The AND-state is readily constructed from 
the Toffoli gate as illustrated in Figure |9j This allows some interesting states to be created 
experimentally, for instance, post-selection of the output to |0) would yield the state 1 00) + 
1 1 ) + 1 10} . (See the course notes [l3| for more on how these techniques can be used as an 
experimental prescription to generate quantum states. See also Figure [10] for the connection 
to computationally universal quantum circuits.) 
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<r 



Quantum 
circuit model 




•/©-units 

|Va> = |000) + |010) + |100) + 1 m> 




FIG. 9. Illustrates the use of units to prepare the AND-state. Using this state together with single 
qubit NOT-gates, one can construct any Boolean qubit state as well as any of the states appearing 
in Table IIV1 We note that the box around the Toffoli gate (left) is meant to illustrate a difference 
between our notation and that of quantum circuits. 



A 



A 





H 









FIG. 10. Hadamard built from the AND-state together with [— } 



= (|0) - |1)). We note that 



def 

quantum computational universality is already possible by considering simple Hadamard states 
(e.g. |^ H ) = |00) + |01) + 1 10> - |11)), COPY- and AND-states, which follows from the proof that 
Hadamard and Toffoli are universal for quantum circuits |67l |. 



1. Summary of the XOR-algebra on tensors 

We will now present the three previously referenced Tables (JH [TT] and IHip which summa- 
rize the quantum logic tensors we introduced in the previous subsections (1III A\ IIIIBI and 
IIIIPp . The tables contain entries listing properties that describe how the introduced net- 
work components interact. These interactions are denned diagrammatically and explained 
in Section HVl 



Gate Type 


Co-copy point (s) 


Unit 


Co-unit Interaction 


COPY 


|0>,|1> 


1+) 


Bell state: |00) + |11> 


Symmetry 


Associative 


Commutative 


Frobenius Algebra 


Full 


Yes 


Yes 


Yes (Node Equivalence) 



TABLE I. Summary of the COPY-gate from Section [III Al 
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Gate Type 


Co-copy point (s) 


Unit 


Co-unit Interaction 


XOR 


!+>,[-) 


|0) 


Bell state: |00) + |11) 


Symmetry 


Associative 


Commutative 


Frobenius Algebra 


Full 


Yes 


Yes 


Yes (Node Equivalence) 



TABLE II. Summary of the XOR-gate from Section [III Bl 



Gate Type 


Co-copy point(s) 


Unit 


Co-unit Interaction 


AND 


|1> 


|1> 


Product state: 11} 


Symmetry 


Associative 


Commutative 


Bialgebra Law 


Inputs 


Yes 


Yes 


Yes (with GHZ) 



TABLE III. Summary of the AND-gate from Section [HI Dl 



E. co-COPY: the co-diagonal 



What is evident from our subsequent discussions on logic gates is that in the context of 
tensors, the bending of wires implies that gates can be used both forwards in backwards. 
We can therefore form tensor networks from Boolean gates in a very different way from 
classical circuits. Indeed, it becomes possible to flip a COPY operation upside down, that is, 
instead of having a single leg split into two legs, have two legs merge into one. In terms of 
tensor networks, co-COPY is simply thought of as being a dual (transpose) to the familiar 
COPY operation. This is common in algebra: to consider the dual notation to algebra, that 
is co-algebra. In general, while a product is a joining or pairing (e.g. taking two vectors and 
producing a third) a co-product is a co-pairing taking a single vector in the space A and 
producing a vector in the space A® A. 



Remark 5 (co-algebras [68]). co-algebras are structures that are dual (in the sense of re- 
versing arrows) to unital associative algebras such as COPY and AND the axioms of which 
we formulated in terms of picture calculi (Sections \IIIA\ and \IIlD\) . Every co-algebra, by 



(vector space) duality, gives rise to an algebra, and in finite dimensions, this duality goes in 
both directions. 

Co-COPY can be thought of as applying a delta function in the transition from input to 
output. That is, given a copy point x — 0, 1, d — 1 for qudits on dim d. Depicting COPY 
as the map A 

A (\x)) = \x) <g> \x) (10) 
we define co-COPY by the map y such that 

V(I»%|J» = W (11) 

that is, the diagram is mapped to zero (or empty) if the inputs \j) do not agree. This 
is succinctly expressed in terms of a delta-function dependent on inputs \i), \j) where i,j = 
0, 1, d — 1 for qudits of dim d. 

Example 6 (Simple co-pairing). Measurement effects on tripartite quantum systems can be 
thought of as co-products. This is given as a map from one system (measuring the first) into 
two systems (the effect this has on the other two). GHZ-states are prototypical examples of 
co-pairings. In this case, the measurement outcome of |0) (|1)) on a single subsystem sends 
the other qubits to 1 00) (|11)) and by linearity this sends |+) to 1 00) + |11). 
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F. The remaining Boolean tensors: NAND-states etc. 

We have represented a logical system on tensors — this enables us to represent any 
Boolean function as a connected network of tensors and hence any Boolean state. We chose 
as our generators, constant |1), COPY, XOR, AND. Other generators could have also been 
chosen such as NAND-tensors. Our choice however, was made as a matter of convenience. If 
we had considered other generators, we could have ended up considering the following cases: 
weak-units (Definition [7]) and fixed point pairs (Definition [9]). We note that the NAND-states 



were used in [69| for fault-tolerant quantum computation — see also [70J . 

Definition 7 (Weak units). An algebra (or product see Appendix [A]) on a tripartite state 
\tp) has a unit (equivalently, one has that the state is unital) if there exists an effect (<f>\ 
which the product acts on to produce an invertible map B, where B = 1 (see Example [H]). If 
no such ((f)\ exists to make B = 1, and B has an inverse, we call (0| a weak unit, and say the 
state is weak unital and if B ^ 1 and B 2 = 1 we call the algebra on \ip) unital-involutive. 
This scenario is given diagrammatically as: 



< 



Example 8 (NAND and NOR). NAND and NOR have weak units, respectively given by |1) 
and |0). These weak units are unital-involutive. 

I^nand) = |001) + |011> + |101> + |110) (12) 

I^nor) = |001) + |010> + 1 100) + |110) (13) 
For IV'nand) to have a unit, there must exist a \<p) such that 

(0|O)|O1) + (0|O)|11) + (0|1)|O1) + (0|1)|1O) = |00) + |11) (14) 

and hence no choice of \<f>) makes this possible, thereby confirming the claim. 

Definition 9 (Fixed Point Pair). An algebra (see Appendix [A]) on a tripartite state 
has a fixed point if there exists an effect (<f>\ (the fixed point) which the product acts on to 
produce a constant output, independent of the other input value. For instance, in Figure 
[T2Tc) on the left hand side the effect (1| induces a map (read bottom to top) that sends 
|+) i — y |1). Up to a scalar, this map expands linearly sending both basis effects (0|, (1| to to 
the constant state |1). If the resulting output is the same as the fixed point, we say (<j)\ has 
a zero (|1) is the zero for the OR-gate in Figure [T2T c)). A fixed point pair consists of two 
algebras with fixed points, such that the fixed point of one algebra is the unit of the other, 
and vise versa (see Figure [TT1 and [T2|) . Diagrammatically this is given in Figure [TT1 



G. Summarizing: network composition of quantum logic tensors 

We have considered sets of universal classical structures in our tensor network model. 
In classical computer science, a universal set of gates is able to express any n-bit Boolean 
function 

/ : B n -> B :: (x u x n ) ^ f( Xl , x n ) (15) 
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< 3- = -^<H «3>_=_ M _ 

FIG. 11. Diagrammatic equations satisfied by a fixed point pair (see Definition [9j) . 




FIG. 12. AND and OR tensors form a fixed point pair. The unit for AND (|1) see a) is the zero for 
OR (c) and vise versa: the unit of OR ([0) see a) is the zero for AND (b). 

where we note that Z 2 = B allowing us to use the alternative notation for / as / : — > 7L d 
with d = 2 for the binary case. Universal sets include {COPY, NAND}, {COPY, AND, NOT}, 
{COPY, AND, XOR,|l)}, {OR, XNOR,|l)} and others. One can also consider the states \ip) 
formed by the bit patterns of these functions f(xi, x 2 ) as 

hM= £ \x 1 )\x 2 )\f(x 1 ,x 2 )) (16) 

a;i,a;2£{0,l} 

This allows a wide class of states to be constructed effectively. In the following Table (llVj) 
we illustrate the quantum states representing the classical function of two-inputs. 



non-linear 


linear (Frobenius Algebras) 


IV'and) = |000) + |010) + |100) + [111) 
I^or) = |001) + |011) + |101) + |111) 
I^nand) = |001) + [011} + |101) + [110} 
I^nor) = [001) + [010} + [100} + [110} 


IVW = |000) + |011) + |101) + |110) 
|V>xnor) = |001) + |010) + |100) + [111} 



TABLE IV. The bit pattern of these quantum states represents a Boolean function (given by the 
subscript) such that the right most bit is the Boolean functions output, and the two left bits are 
the functions inputs, and the nondinear Boolean functions are on the left side of the table and the 
linear functions on the right. Consider the state IV'and)! an d Boolean variables x\ and x 2 , then the 
superposition (Vand) encodes the function \xi,x 2 ,x\ A x 2 ) in each term in the superposition, and 
IV'and) = Yl Xl z 2 e{o 1} \ x li x 2> x i A x 2 ). As outlined in the text, cup/cap induced-duality allows us 
(for instance) to express this state as the operator |0)(00[ + 1 0) (01 [ + |0)(01| + |1)(H| :: \x\,x 2 ) 
\x\ A x 2 ) which projects qubit states to the AND of their bit value. 
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IV. INTERACTION OF THE NETWORK COMPONENTS 

Having outlined the Boolean components used in our tensor toolbox we now explore 
how these tensors interact when connected in a tensor network. The interactions can be 
defined diagrammatically and given simple rewrite rules for CTNS based on these component 
tensors. 



1. Merging COPY-dots by node equivalence 



COPY-dots are readily generalized to an arbitrary number of input and output legs. As 
one would rightly suspect, a COPY-dot with n inputs and m outputs corresponds to an 
n + m-partite GHZ-state. Neighboring dots of the same type can be merged into a single 
dot: this is called node equivalence in digital circuits. COPY-dots represent Frobenius 



algebras [68|, |7l|, l22| . (Note that the recent online version of 



the arXiv version of the present manuscript as well as [13j, [66 ] . 



72[ was already influenced by 



Theorem 10 (Node equivalence or spider law). Given a connected graph with m inputs and 
n outputs comprised solely of COPY dots of equal dimension, this map can be equivalently 
expressed as a single m-to-n dot, as shown in Figure VTl 



Example 11 (Two-site reduced density operator of n-party GHZ-states). GHZ-states on 
n-parties have a well known matrix product expression given as 



|GHZ n ) = Tr 



/|0) 
|1> 



LV o 



o 





\d-l)J J 



|0)|0)|0) + |l)|l)|l) + ... + |d-l)|d-l)|d-l) 



(17) 

Such MPS networks are known to be efficiently contactable. We note that the networks in 
Figure [13] do not appear a priori to be contractible due to the number of open legs. What 
makes them contractible (in their present from) is that the tensors obey node equivalence law 
allowing them to be deformed into a contractible MPS network: see Figure HU The reduced 
density matrix of an n-party GHZ-state then becomes (a) in Figure [15] and the expectation 
value of an observable is shown in (b) where we included the normalisation constant. 




FIG. 13. Node equivalence or spider law [68]. Connected black-dots (•) as well as connected plus- 
dots (©) can be merged and also split apart at will. The intuition for digital or qudit circuits follows 
by connecting a state \<f>) to one of the legs and iterating over a complete basis |0), — 1). 
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FIG. 14. The GHZ-state tensor is simply a rank-ra COPY-dot. Node equivalence implies that this 
tensor can be deformed into any network geometry including a MPS comb-like structure (right). 




= Tr(p' GHZ -O a ®0 2 ) 



FIG. 15. Reduced density operator. Left (a) reduced density operator p' GHZ found from applying 
the node equivalence law to a n-qubit GHZ-state. Right (b) the expectation value of observable 
0\ ®0 2 found from connecting the observable and connecting the open legs (i.e. taking the trace). 



A. Associativity, distributivity and commutativity 

The products we have considered are all associative and commutative. As algebras, AND, 
XOR and COPY are associative, unital commutative algebras. This was already expressed 
diagrammatically in Figures [3]^a) and Figure H^c). The diagrammatic laws relevant for this 
subsection represent the following Equations 

(xi A x 2 ) A x 3 = xi A (x 2 A x 3 ) (18) 

(xi © x 2 ) © x 3 = Xi © (x 2 © x 3 ) (19) 
Distributivity of AND over XOR then becomes (see (h) in Figure H]) 

(xt © x 2 ) Ax 3 = (xi A x 2 ) © (xi A x 2 ) (20) 

We have commutativity for any product symmetric in its inputs: this is the case for AND 
and XOR. 



B. Bialgebras on tensors 

There is a powerful type of algebra that arises in our set ting : a bialgebra denned graph- 
ically on tensors in Figure US (see Kassel, Chapter III @, H). 



Such an algebra is simultaneously a unital associative algebra and co-algebra (for the 
associativity condition see (b) in Figure [TB]) . Specifically, we consider the following two 
ingredients: 

(i) : a product (black dot) with a unit (black triangle) see the right hand side of Figure [ToT a) 

(ii) : a co-product (white dot) with a co-unit (white triangle) see the left hand side of 

Figure [TBTa) 



To form a bialgebra, these two ingredients above must be characterized by the following four 
compatibility conditions: 
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(i) : The unit of the black dot is a copy-point of the white dot as in (e) from Figure [T6l 

(ii) : The (co)unit of the white dot is a copy-point of the black dot as in (d) from Figure fT6l 

(iii) : The bialgebra-law is satisfied given in (c) from Figure [161 

(iv) : The inner product of the unit (black triangle) and the co-unit (white triangle) is 

non-zero (not shown in Figure [T6|) . 



w ti_| a ,bl ^^ ,C, U U <d, M m «>$ 



FIG. 16. Bialgebra axioms 68| (scalars are omitted), (a) unit laws (these are of course left and 



right units); (b) associativity; (c) bialgebra; (d,e) co-COPY points. 

Example 12 (GHZ, AND form a bialgebra). We are in a position to study the interaction 
of GHZ-AND. This interaction satisfies the equations in Figure [TBI (a) the bialgebra law; (b) 
the co-copy point of AND is |1); and (c) the co-interaction with the unit for GHZ creates a 
compact structure. In addition, (a) and (b) show the copy points for the black GHZ-dot; in 
(c) we have the unit and fixed point laws. 

Even if a given product and co-product do not satisfy all of the compatibility conditions 
(given in (a), (b), (c), (d), (e) in Figure fTBT). and hence do not form bialgebras, they can 
still satisfy the bialgebra law which is given in Figure [ToT c). Examples of states that satisfy 
the bialgebra law in Figure IToT c). but are not bialgebras are given in Definition [T31 Notice 
that bialgebra provides a highly constraining characterization of the tensors involved and is 
tantamount to defining a commutation relation between them. 



Definition 13 (Bialgebra Law [681]). A pair of quantum states (black, white dots) satisfy 
the bialgebra law if (c) in Figure \M holds. The Boolean states, AND, OR, XOR, XNOR, 
NAND, NOR all satisfy the bialgebra law with COPY. 



1. Hopf algebras on tensors 



A particularly important class of bialgebras are known as Hopf-algebras [68J. This is 
characterized by the way in which algebras and co- algebras can interact. This is captured 
by the Hopf-law, where the linear map A is known as the antipode. 

Definition 14 (Hopf-Law [68]). A pair of quantum states satisfy the Hopf-Law if an A can 
be found such that the following equations hold: 




wjspf-Law 




Example 15 (XOR and COPY are Hopf-algebras on Boolean States [35|] ) . It is well known 
(see e.g. [35[) that the Boolean state XOR, satisfies the Hopf-algebra law with trivial antipode 
(A = 1) with COPY. Recall Figure E2(g) . 
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C. Bending wires: compact structures 



As mentioned in the preliminary section ((IT]), we make use of what's called a compact 
structure in category theory which amounts to introducing cups and caps, to provide a 
formal way to bend wires and define transposition. See Figures [17] and [18] 

A compact structure on an object % consists of another object H* together with a pair 
of morphisms (note that we use the equation %* = % in Hilbert space making objects self 
dual which simplifies what follows). The theory is well known in the modern mathematics 



of category theory and has been used to study teleportation in [74|]. Consider 



m ■ 1 



e H :H®H 



where the standard representation in Hilbert space with dimension d and basis {|i}} is given 
by 

d-l d-l 

VH = y~"; \i) ® \i) e H = ^ (i| ® (i\ 

8=0 1=0 

and in string diagrams (read from the top to the bottom of the page) as 



(a) 



(b) 



These cups and caps give rise to cup/cap-induced duality: this amounts to being able to 
create a linear map that "flips" a bra to a ket (and vise versa) and at the same time taking 
an (anti-linear) complex conjugate. In other words, the cap J2l=o sen ds quantum state 
= a\0) + j3\l) to a(0| + f3(l\ which is equal to the complex conjugate of \ip)^ = = 
a(0\ + /3(1\. Diagrammatically, the dagger is given by mirroring operators across the page, 
whereas transposition is given by bending wire(s). Clearly, {ip\ = a(0\ + 

In the case of relating the Bell-states and effects to the identity operator, under cup/cap- 
induced duality, we flip the second ket on r] H and the first bra on e%. This relates these 
maps and the identity 1% of the Hilbert space: that is, we can fix a basis and construct 
invertible maps sending r/y^ = 1% = e%. More generally, the maps 77% and e% satisfy the 
equations given in Figure [17] and their duals under the dagger. 

A second way to introduce cups and caps is to consider a Frobenius form 



68 on either 



of the structures in the linear fragment from Figure [3] (COPY and XOR). This is simply a 
functional that turns a product/co-product into a cup/cap. This allows one to recover the 
above compact structures (that is, the cups and caps given above) as 



(a) 




(b) 



Y 



Again, we will use these cups and caps as a formal way to bend wires in tensor networks: 
this can be thought of simply as a reshape of a matrix. 
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Ca) 




(di 



c 



f - 




FIG. 17. Cup identities, (a) Symmetry, (b) Conjugate state, (c) Teleportation 6.1] or the snake 
equation, (d) Sliding an operator around a cup transposes it. 




(b) 



FIG. 18. Diagrammatic adjoints. Cups and caps allow us to take the transpose of a linear map. 
Note that care must be taken, as flipping a ket \tp) to a bra (ip\ is conjugate transpose, and bending 
a wire is simply transposition, so the conjugate must be taken: e.g. acting on \ip) with a cap given 
as Y2i {m\ results in 



V. EXAMPLES OF CATEGORICAL TENSOR NETWORK STATES 

Our categorical approach enables one to translate a quantum state directly into a new 
type of network: a so-called CTNS. We have focused on Boolean network components and 
have already presented in detail their algebraic properties and defining characteristics. Here 
we will illustrate their expressive power by considering a few elementary examples before 
presenting our main theorem (1231) . precisely showing how to determine a categorical tensor 
network to represent any given quantum state. 



A. Constructing Boolean states 

Since the fixed building blocks of our tensor networks are the logic tensors AND, OR, 
XOR and COPY, along with ancilla bits, we can immediately apply the universality of 
these elements for classical circuit construction to guarantee that any Boolean state has 
a categorical tensor network decomposition. However our construction goes beyond this 
because as we have seen, categorical tensor networks can be deformed and rewired in ways 
which are not ordinarily permitted in the standard acyclic-temporal definition of classical 
circuits. The W-state will be shown to provide a non-trivial example of this. 

Example 16 (Functions on W- and GHZ-states). We consider the function /w which outputs 
logical-one given input bit string 001, 010 and 100 and logical-zero otherwise. Likewise the 
function f GHZ is defined to output logical-one on input bit strings 000 and 111 and logical- 
zero otherwise. See Examples [T^l and I2"0l which consider representation of these functions as 
polynomials. We will continue to work with a linear representation of functions on quantum 
states; here bit string 000 1 000) (etc.). 
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Example 17 (MPS form for W-state). Like the GHZ state, the W-state has a simple MPS 
representation 

l W -> = <°l ( !?\ ln\ 1 I 1 ) = l 10 -°) + l 01 -°) + - + I 00 - 1 )- (21) 



.11) |0>. 

This description (I2T1) is succinct. All MPS-states have essentially this same topological or 
network structure. In contrast, our categorical construction described below breaks this 
network up further. 

Remark 18 (Exact- value functions). The function /w takes value logical-one on input vec- 
tors with k ones for a fixed integer k. Such functions are known in the literature as Exact- 
value symmetric Boolean functions. When cast into our framework, exact-value functions 
give rise to tensor networks which represent what are known as Dicke states fSaj . 

Example 19 (Function realization of /w and /ghz: the Boolean case). One can express 
(using x to mean Boolean variable negation) 

fw(xi, x 2 , x 3 ) = x{x 2 xz © x{x 2 x 3 © xix 2 x 3 (22) 

by noting that each term in the disjunctive normal form of /w are disjoint, and hence 
OR maps to XOR as V i— > ©. The algebraic normal form (see Appendix [B]) becomes 

f\N{x\,X 2 , X 3 ) = Xi © X 2 © X 3 © XiX 2 X 3 (23) 
/GHZ(^1, %2, X 3 ) = 1 © Xi © X 2 © X 3 © XiX 2 © XiX 3 © x 2 x 3 (24) 

Example 20 (Function realization of /w and /ghz: the set function case). Set functions 
are mappings from the family of subsets of a finite ground set (e.g. Booleans) to the real or 
complex numbers. In the circuit theory literature, functions from the Booleans to the reals 
are known as pseudo-Boolean functions and more commonly as multi-linear polynomials or 
forms (see [75] where these functions are used to embed a co-algebraic theory of logic gates in 
the ground state energy configuration of spin models). There exists an algebraic normal form 
and hence a unique multi-linear polynomial representation for each pseudo-Boolean function 
(see Appendix [B]) . This is found by mapping the negated Boolean variable asi4 (1 —x). 
For the GHZ- and W-functions defined in Example [16] we arrive at the unique polynomials 
025]) and 026]). 

/ghz(#i, x 2) x 3 ) = 1 — xi — x 2 + x\x 2 - x 3 + xtx 3 + x 2 x 3 (25) 

fw(xi, x 2 , x 3 ) = xx + x 2 + x 3 - 2x\x 2 - 2x\x 3 - 2x 2 x 3 + 3xix 2 x 3 (26) 

These polynomials (1231) and ( 1261) are readily translated into categorical tensor networks. 

Example 21 (Network realisation of W- and GHZ-states). A network realization of W- and 
GHZ-states in our framework then follows by post-selecting the relevant network to |1) on 
the output bit — leaving the input qubits to represent a W- or GHZ-state respectively. An 
example of this is shown in Figure [Tj|] 

Two different categorical constructions for the building blocks of the W-state are shown 
in Figure [2D] and Figure [2TJ Notice that in Figure [2T] the resulting tensor network forms 
an atemporal classical circuit and is much more efficient than the naive construction in 
Figure [2D] Moreover by appropriately daisy-chaining the networks in Figure [2T] we construct 
a categorical tensor network for an n-party W-state as shown in Figure [22] The resulting 
form of this tensor network is entirely equivalent (up to regauging) to the MPS description 
given earlier, but now reveals internal structure of the state in terms of CTNS building 
blocks. 
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fw(x 1 ,x 2 ,x 3 ) 



(b) 
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FIG. 19. Left (a) the circuit realization (internal to the triangle) of the function /w of e.g. f|23|) which 
outputs logical-one given input \x1X2X3) = |001), 1 010) and 1 100) and logical-zero otherwise. Right 
(b) reversing time and setting the output to |1) (e.g. post-selection) gives a network representing 
the W-state. The naive realization of /w is given in Figure [21] with an optimized co-algebraic 
construction shown in Figure [2TJ 




FIG. 20. Na'ive CTNS realization of the familiar W-state |001) + |010) + |100). A standard 
(temporal) acyclic classical circuit decomposition in terms of the XOR-algebra realizes the function 
/w of three bits. This function is given a representation on tensors. As illustrated, the networks 
input is post selected to |1) to realize the desired W-state. 



B. Describing states with complex coefficients 

Boolean states, such as the GHZ- and W-states, are typified by being superpositions 
of computational basis states with equal real coefficients (in both cases, these coefficients 
take only binary values, and 1). In this section, we will permit a minor extension to 
binary superposition input/output states by considering arbitrary rank-1 tensors within our 
otherwise Boolean tensor networks. This is illustrated by a simple example: 

Example 22 (Network realization of = |01) + 1 10} + cvjt|ll)). We will now design a 
network to realize the state |01) + 1 10} + «fc|ll). The first step is to write down a function 
fs such that 

/ S (0,l) = / 5 (1,0) = / 5 (1,1) = 1 (27) 

and /s(00) = (in the present case, fs is the logical OR-gate). We post select the network 
output on |1), which yields the state 1 1 } + 1 10} + |11), see Figure [231(a). The next step is to 
realize a diagonal operator, that acts as identity on all inputs, except 1 1 1 } which gets sent 
to a*; 1 11). To do this, we design a function f d such that 



/ d (0,l) = / d (l,0) = / d (0,0) = 



(28) 
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110) + 1 101) + |011) 



FIG. 21. W-class states in the categorical tensor network state formalism, (a) is the standard 
W-state. (b) is found from applying De Morgan's law (see Section IIII Dj) to (a) and rearranging 
after inserting inverters on the output legs. Notice the atemporal nature of the circuits, as one 
gate is used forwards, and the other backwards. 




FIG. 22. W-state (n-party) in the categorical tensor network state formalism. The comb-like feature 
of efficient network contraction remains, with the internal structure of the network components 
exposed in terms of well understood algebraic structures. 

and fd(l, 1) = 1 (in the present case, fd is the logical AND-gate). This diagonal, takes the 
form in Figure 123T b). The final state \ip) = |01) + 1 10) + afc|ll) is realized by connecting 
both networks, leading to Figure 1237 c). 




oi> + io> + 



FIG. 23. Categorical tensor network representing state = |01) + 1 10) + afc|ll), as explained in 
Example El 



VI. PROOF OF THE MAIN THEOREMS 

We are now in a position to state the main theorem of this work. Specifically, we have a 
constructive method to realize any quantum state in terms of a categorical tensor network. 
(A corollary of our exhaustive factorization of quantum states into tensor networks is a new 
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type of quantum network universality proof. To avoid confusion, we point out that past uni- 
versality proofs in the gate model already imply that the linear fragment (Figure |3]) together 
with local gates is quantum universal. However, the known universality results clearly do 
not provide a method to factor a state into a tensor network! Indeed, the decomposition 
or factorization of a state into a tensor network is an entirely different problem which we 
address here.) We state and prove the theorem for the case of qubit. The higher dimen- 
sional case of qudits follows from known results that any d-state switching function can be 



expressed as a polynomial and realized as a connected network [43|, |76|, [77J . The theorem 
can be stated as 

Theorem 23 (Tensor network representation of quantum states). For any state \ip) of 
n-qubits with the form 

k 

lV>) = 5>;l0i), (29) 

where otj are complex coefficients and for each j the state \<pj) is an equal superposition of 
a set of computational basis states, it can be represented as a network containing tensors 
from the quantum Boolean calculus (Figures [3] anrfJ^P, together with input/ ouput states of 
the form \aj) := |0) + Oj|l). 

Notice that an arbitrary state can be brought into the form required of \ip) by composing 
it as k = 2 n terms with each state \<f>j) being a single distinct computational basis state. 
The proof is simplified by invoking some supporting lemmas. 

Lemma 24. There exists a map g represented by a tensor network taking diagonal maps in 
rt C 2 -> <g) n C 2 onto quantum states in ® n C 2 . 

Proof. Let V be a diagonal map in (^) n C 2 — >■ n C 2 . We write V = 2^ x e{o iv» a x|x)(x| and 
proceed as follows (where the term £>°(S) n (|0) + immediately yields the desirable tensor 
network depiction) 



g{V}:=Vo0(\O) + \l))=Vo \V)= E «x|x)(x|y) = 

n ye{o,i} n x,ye{o,i}™ 

= E <5xy«x|x) = X a x |x) 
xe{o,i}™ xe{o,i}" 

□ 

Lemma 25. There exists a map h represented by a tensor network taking quantum states 
in (§Z) n C 2 onto diagonal maps in C 2 — > (^) n C 2 . 

Proof. Let \ip) be a quantum state in (^) n C 2 . Let V be a diagonal map in (^) n C 2 — > n C 2 . 
We write V = J2 x e{o i} n a x|x)(x| and proceed as follows (where the term (§) n (X/i=o/i 
Sxe{o i}™ a x|x) immediately yields the desirable tensor network depiction) 

fcW»:=<8>(£ E a ^= E E a *ix> = 

n i=o/i xe{o,i}" ye{o,i} n *.e{o,i} n 

= E lyy) a x^xy = E a *l xx ) 

x,y6{0,l}™ x£{0,l}" 
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and then we now write h in terms of h! 

MIV'>}:=®(E< << I)° /I, w = E (yyi° E a *i xx > 

n i=o/i ye{o,i}™ xe{o,i}™ 

= ^2 a x|x)(y|5 yx = ^2 «x|x)(x| 

y, xe{o,i}™ xe{o,i}™ 



□ 



Corollary 26. It follows that g{h{\^)}} — l^o \tp) = and h{g{V}} = lpoD = V and 

hence we have inverses for g and h establishing an isomorphism between diagonal operators 
in n C 2 — > (£) n C 2 and quantum states in (£) n C 2 . 

With the supporting lemmas in place, we will now proceed to prove Theorem l23l 

Proof. Returning to our particular expression for an arbitrary quantum state (note that in 
Equation ( I50|) we now append an extra term, by letting the sum run from j — instead of 
j — 1, the use of this will become clear below) 

k 

iv>) = E a ^>> ( 3 °) 

3=0 

where ctj are complex coefficients and for each j the state \<f>j) is an equal superposition of 



arc 



a set of computational basis states, we will explain how k + 1 asynchronous circuits [41 
used to factor the state, and express the state as a CTNS (here and in what follows k is the 
highest term in the sum from Equation (j29p ). 

We proceed by returning to our original expression (1291) from Theorem [23] (starting from 
j = 1) with the coefficients removed 

w) = E ifc> ( 31 ) 

3=1 

Each individual term in Equation fl3Tl) is then expressed in the computational basis and 
used to form a set denoted C + . All corresponding bit patterns of the same dimension 
not appearing in this expression form a second set C~ (where clearly C + fl £~ = and 
C + U C~ = {0, 1}™, where n is the number of qubits in the desired state). We proceed to 
construct a function f that outputs logical-one on all input bit strings in C + and outputs 
local zero on all input bit strings in L~ . The function acts on n + 1 bits, the inputs are given 
on the right of the tensor symbol and the output on the left of the tensor symbol in (|32|) 

\4>f)= E l x )®l/W) ( 32 ) 
xe{o,i} n 

where /(x) : B n — > B :: x i— > /(x) was given the representation on quantum states in Section 
IVl Post selecting the networks output (the rightmost bit in Equation f[52"j) ) |1) realizes the 
desired superposition of terms, with all coefficients of the terms and hence relative phases 
equal. 

W f )= E WWW (33) 
xe{o,i} n 
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For our specific construction, depicted in Figure [2U we proceed by inverting the output 
of the function (e.g. f h-> f Q © 1). We then post select the output of the function to the 
state \a ) = |0) + « |1) = |0) for the choice a = 0. This handles the j = term in the sum 
Equation (130]) . 

To adjust the amplitudes of the desired state from Equation ( 1301) . we will construct 
tensors that represent diagonal operators. For the jth term in with coefficient ctj, we 
again construct a function fj. We represent \<f>j) in the computational basis, and each term 
in this expression is used to form a set denoted C + . All corresponding bit patterns of the 
same dimension not appearing in this expression form a second set C~ . We proceed to 
construct fj to output logical-one on all input bit strings in C + and outputs logical-zero on 
all input bit strings in Cr . The network is then post selected to |0) + which results in 
states of the form 

hM= £ (0|/(x))|x) + a j J2 WW)!*) ( 34 ) 
xe{o,i}" xe{o,i} n 

and we transform fj into a diagonal operator having entries G {1, a,} by applying the map 
h from Lemma [251 resulting in the diagonal map 

V 3 = (0|/(x))|x)(x| + % (m*))\x)(x\ (35) 

xe{0,l} n x6{0,l} n 

We will apply k such commuting maps T>j to the initial state, accounting for k asynchronous 
circuits. The operators are composed by means of n co-COPY-dots from Section Till El (see 
Figure Upland Example |2"2"|) . There will be a single output with open legs which gives the 
state. Each of the n COPY-dots will then require k + 2 legs. The resulting construction then 
gives tensor networks with the form shown in Figure [2H □ 




FIG. 24. The CTNS for a state resulting from our exhaustive construction procedure. 
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Remark 27 (Qudit states). In TheoremW^ we considered the arbitrary states ofn-qubits. By 
using multi-valued logic (also called d-state switching, or many-valued logic), it is possible to 
define a universal gate set similar to what was done for the case of qubits and so equivalently 
construct a CTNS for n-body qudit systems f^l . Xl& , VF\ l. 

Definition 28 (Generalized Polynomial Boolean States (GPBS)). The size of a tensor 
network is the number of tensors it contains and its depth is the maximal length of a 
path from any tensor to any other. Consider families of uniform circuits, built from the 
XOR-algebra, that is the bounded fan-in gates AND, XOR, COPY, of arity two and the 
constant |1). We will index these circuit families by bounding the circuit depth, which 
also has the impact of bounding the maximum fan-in and the circuit size. We will then 
consider categorical tensor networks to represent states of n interacting d- level systems 
constructed from bounded circuit families to realize each fj with the form given in Figure 
¥2M We will proceed by indexing these families of categorical tensor networks in terms of k, 
the maximum depth of any given circuit realizing any function fj in the network. We will 
then bound the number of such functions fj to be at most some polynomial in k. We then 
determine how k(n) changes. This works by considering circuit families and categorical 
tensor networks of a given form, used to represent quantum states on increasingly many 
subsystems n. If k(n) is bounded by a polynomial in n the categorical tensor network has 
an efficient description. We index such families as C(/c), and refer to them as Generalized 
Polynomial Boolean States (GPBS). 

Theorem 29. GPBS from Definition are sampled exactly in the computational basis in 
time and space complexity bounded by poly(fc). 

Proof. To prove Theorem [291 we begin first by considering a qudit state vector \xq, Xi, x n ) 
for specific xq, xi, x n G {0, 1, d — 1}. We wish to know the coefficient (xo, x%, x n \C), 
where C is a CTNS representing a GPBS. The COPY-gates in the construction from Figure 
[2^1 map 

\x ,Xi, ...,x n ) >->• \x ,x 1 , ...,£„) (36) 

poly(fe) 

each of these poly(A;) vectors will be acted on by a network realizing fj post selected to the 
state \a.j). Hence, to sample the network C amounts evaluating the sum 

(x ,xi,...,x n \C)= (fj(x ,xi, ...,x n )\aj) = c E C (37) 

jSpoly(fc) 

The proof then follows by simply evaluating each of the poly(fc) poly(/c)-depth functions 
fj(xo,Xi, ...,x n ) and then summing the inner products (fj(xo,Xi, ...,x n )\aij). □ 

We note that the construction in Theorem [23] automatically groups basis states with the 
same coefficients ctj, of the k terms. Further reductions are also possible if say a given set of 
coefficients are given by products of other coefficients. While this construction does prove 
the existence of a CTNS (along with how to build it) our construction will not render efficient 
representations for general cases, as one might expect. Indeed, there is no guarantee that 
any of the k + 1 switching functions are efficient in their complexity, nor that the resulting 
complete network is contractible. The latter property is in fact confounded by the presence 
of fan-in (up to fc + 2 legs) of the n co-COPY-dots (the presence also implies that the network 
cannot represent a deterministic physical preparation procedure [66|). However, as we saw 
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earlier with string-bond states in Figure [2j the COPY dot breaks up when computational 
basis states are inputted. For our general decomposition in in Figure EU this case causes 
the k + 1 Boolean switching functions to similarly decouple. Intuitively, if we further restrict 
ourselves to k + 1 being polynomial in n, and additionally that each switching function has 
a depth which is also polynomial in n, then the subsequent evaluation of the amplitude of 
the state is efficient for any computational basis state (see Definition [28] and Theorem 1291) . 
This is a weak requirement in practice and interestingly, does not depend on the internal 
geometry of the networks representing the functions, but only on their depth and size. Thus 
we have found a new general class of states which can be sampled exactly and efficiently. 
Finally the construction was based on using acyclic-temporal Boolean circuits. However, we 
have already seen that in the tensor context wires can be bent around: it is not necessary 
for a tensor network to correspond to a valid classical circuit. As the W-state example 
(see Figures EH and I2"2"j) illustrated the tensor networks can be much simpler once this new 
freedom is exploited. 



VII. OUTLOOK AND CONCLUDING REMARKS 



We have introduced methods from category theory and algebra to tensor network states. 
Our main focus has been on tensor networks built from contracting Boolean tensors and 
we have outlined their algebraic properties. The logical conclusion of this approach has 
led us to a exhaustive CTNS decomposition of quantum states. From this we obtained 
a class of quantum states, which can be sampled in the computational basis efficiently 
and exactly. The expressiveness and power of this new method was further illustrated by 
considering several simple test cases: we considered internal structure of some MPS states, 
e.g. GHZ- and W-states. We have opened up some future potential research directions. In 
particular, beyond the form of our construction there is an open question as to whether the 
algebraic properties of some subset of the tensors in CTNS can enable efficiently contractible 
networks beyond those already known which are based on topology (like MPS) and additional 
unitary/isometric constraints (MERA). In this way future studies of CTNS may lead to new 
classes of states and algorithms which will help challenge and shape our understanding of 
many-body physics. 
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Appendix A: Algebra on quantum states 



We are concerned with a network theory of quantum states. This on the one hand can 
be used as a tool to solve problems about states and operators in quantum theory, but does 
have a physical interpretation on the other. This is not foundational per se but instead 
largely based on what one might call an operational interpretation of quantum states and 
processes. We call an algebra a pairing on a vector space, taking two vectors and producing 
a third (you might instead call it a monoid if there is a unit, and then a group if the set 
of considered vectors is closed under the product). Let's now examine how every tripartite 
quantum state forms an algebra 13 . 

Consider a tripartite quantum state (subsystems labeled 1,2 and 3), and then ask the 
question: "how would the state of the third system change after measurement of systems 
one and two?" Enter Algebras: as stated, an algebra on a vector space, or on a Hilbert 
space is formed by a product taking two elements from the vector space to produce a third 
element in the vector space. Algebra on states can then be studied by considering duality 
of the state, that is considering the adjunction between the maps of type 



1 ^U®U®U 



and 



(Al) 



This duality is made evident by using the f-compact structure of the category (e.g. the cups 
and caps). It is given vivid physical meaning by considering the effect measuring (that is 
two events) two components of a state has on the third component. 

Remark 30 (Overbar notation on Spaces). Given a Hilbert space H, we can consider the 
Hilbert space % which can be simply thought of as the Hilbert space ~H will all basis vectors 
complex conjugates (overbar). That is, H, is a vector space whose elements are in one-to-one 
correspondence with the elements ofH: 



H = {v | v e %}, 

with the following rules for addition and scalar multiplication: 



v + w = v + w and av = av 



(A2) 



(A3) 



Remark 31 (Definition of Algebra). We consider an algebra as a vector space A endowed 
with a product, taking a pair of elements (e.g. from A® A) and producing an element in A. 
So the product is a map A <g> A — > A, which may not be associative or have a unit (that is, 
a multiplicative identity — see Example^ for an example of an algebra on a quantum state 
without a unit). 

Observation 32 (Every tripartite Quantum State Forms an Algebra). Let \ip) G H <S> H <8> 
H be a quantum state and let M i; Mj be complete sets of measurement operators. Then 
forms an algebra. 



GV'), Mi, Mj 




^/-product 




= 4- 
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The quantum state = ^2 i j k '4 ) ^ k \ijk) is drawn as a triangle, with the identity operator 
on each subsystem acting as time goes to the right on the page (represented as a wire). 
Projective measurements with respect to M, and Mj are made. We define these complete 
measurement operators as 

TV 

M 1 = ^2i\tf> i )(i/} i \ (A4) 
i=i 



N 

M 2 = Y,j\h)(H (A5) 

3=1 

such that we recover the identity operator on the A/"- level subsystem viz 

N N 

£l^>^;l = Z>«>M = 1 * (A6) 
j=i i=i 

The measurements result in eigenvalues i,j leaving the state of the unmeasured system in 

\tu) = Y,r vz (T\x)(^\y)\z) (A7) 

xyz 

where (Q\ = \Q) T that is, the transpose is factored into: (i) taking the dagger (diagram- 
matically this mirrors states across the page) and (ii) taking the complex conjugate. Hence, 

= |Q) T = {Q\=W (A8) 
and if we pick a real valued basis for \x), \y), \z) = |0), |1) we recover 

\co) = j2r yz (xm(y\<f>y}\z) (A9) 

xyz 

As stated, this physical interpretation is not our main interest. Even in its absence, we're 
able to write down and represent a quantum state purely in terms of a connected network, 
where each component is fully defined in terms of algebraic laws. 



Appendix B: XOR-algebra 

Here we review the concept of an algebraic normal form (ANF) for Boolean polynomials, 
commonly known as PPRMs. See the reference book 63] and the historical references [64], 65| 
for further details. 

Definition 33. The XOR-algebra forms a commutative ring with presentation M = 
{B, A, ©} where the following product is called XOR 

— :BxB^B:: (a,b) ^ a + b - ab mod 2 (Bl) 
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and conjunction is given as 

- A— : B x B B :: (a,b) -»■ a ■ b, (B2) 

where a ■ b is regular multiplication over the reals. One defines left negation ->( — ) in terms 
of © as -.(— ) = 

1 © ( — ) :Bn-B::a-^l-a. (B3) 

In the XOR-algebra, 1-5 hold, (i) a © = a, (ii) a © 1 = ->a, (iii) a © a = 0, (iv) a ffi -ia = 1 
and (v) a V 6 = a © 6 © (a A 6). Hence, is the unit of XOR and 1 is the unit of AND. The 
5th rule reduces to aV6 = a©6 whenever a A b = 0, which is the case for disjoint (mod 2) 
sums. The truth table for AND follows 



Xi 


X 2 


f(xi,x 2 ) = xi A x 2 














1 





1 








1 


1 


1 



Definition 34. Any Boolean equation may be uniquely expanded to the fixed polarity 
Reed-Muller form as: 

fyXi, X 2 , • Xf.) = Co © CiX^ 1 © C 2 X^ 2 ffi ' " " © C n X n n © 

c n+1 x?x° n © • ■ ■ © c 2k . x x\ x x a 2 \ 4 fe , (B4) 

where selection variable G{ e {0, 1}, literal represents a variable or its negation and any c 
term labeled Cq through Cj is a binary constant or 1. In Equation flB4p only fixed polarity 
variables appear such that each is in either un-complemented or complemented form. 

Let us now consider derivation of the form from Definition |3H Because of the structure 
of the algebra, without loss of generality, one avoids keeping track of indices in the node 
case, by considering the case where A^ = 2 n = 8. 

Example 35. The vector c = (co, c±, c 2 , C3, C4, C5, c$, C7, ) T represents all possible outputs of 
any function f(x\, x 2 , X3) over the algebra formed from linear extension of Z 2 x Z2 x "Z 2 . We 
wish to construct a normal form in terms of the vector c, where each c, G {0, 1}, and therefore 
c is a selection vector that simply represents the output of the function / :BxBxB-> 
B :: (xi,x 2 ,x 3 ) 1— )■ f(xi,x 2 ,x 3 ). One may expand / as: 

f(x!, x 2 , x 3 ) = (co ■ -ixt ■ -.x 2 • ^r 3 ) V (ci • -ia?i ■ • ar 3 ) V (c 2 ■ ->Xi • x 2 ■ ->x 3 ) 
V(c 3 • -1x1 ■ x 2 ■ x 3 ) V (c 4 ■ Xi • -ix 2 • -ix 3 ) V (c 5 • X\ ■ -ix 2 • x 3 ) 
V(c 6 • xi • x 2 • ->x 3 ) V (c 7 • xi • x 2 • x 3 ) (B5) 

Since each disjunctive term is disjoint the logical OR operation may be replaced with 
the logical XOR operation. By making the substitution ->a = a © 1 for all variables and 
rearranging terms one arrives at the following normal form. (For instance, ->Xi • -1X2 • ~^x 3 = 
(1 © x x ) ■ (1 © x 2 ) ■ (1 © x 3 ) = (1 © x 1 © x 2 © x 2 ■ x 3 ) ■ (1 © X 3 ) = 1 © Xi © X 2 © X 3 © Xi ■ x 3 © 
x 2 ■ X 3 © Xi ■ x 2 • x 3 .) 
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f(xi, x 2 , x 3 ) = c © (Co © c 4 ) ■ Xi © (c © c 2 ) ■ X 2 © (Co © Ci) ■ x 3 © (c © c 2 © c 4 © C 6 ) ■ aq • x 2 
ffi(c © Ci © C 4 © C 5 ) • X\ ■ x 3 © (c © Ci © c 2 © c 3 ) • x 2 ■ x 3 
©(c © ci © c 2 © c 3 © c 4 © c 5 © c 6 © oi) -xi-x 2 -x 3 (B6) 

The set of linearly independent vectors, {xi, x 2 ,x 3 , aq ■ x 2 ,X\ ■ x 3 ,x 2 ■ x 3 ,X\ ■ x 2 ■ x 3 } 
combined with a set of scalars from Equation IB6I spans the eight dimensional space of the 
Hypercube representing the Algebra. A similar form holds for arbitrary N. 

f(xi,x 2 , x 3 ) = (ai) • X\ © (a 2 ) • x 2 © (x 3 ) ■ x 3 © (ai © a 2 © ai © c 2 ) ■ Xi ■ x 2 
ffi(ai © a 3 © ai © c 3 ) • xi ■ x 3 © (a 2 © a 3 © a 2 © c 3 ) • x 2 • x 3 
ffi(ai © a 2 © a 3 © ai © a 2 © a 3 ) • • x 2 ■ x 3 (B7) 



[1] Saunders Mac Lane. Categories for the working mathematician 2nd ed. Graduate Texts in 

Mathematics, Springer, 1998. 
[2] S. J. Denny, J. D. Biamonte, D. Jaksch, and S. R. Clark. Algebraically contractible topological 

tensor network states. J. Phys. A: Math. Theor., 45:015309, 2012. arXiv:1108.0888. 
[3] S. Al-Assam, S. R. Clark, C. J. Foot, and D. Jaksch. Capturing long range correlations in 

two-dimensional quantum lattice systems using correlator product states. ArXiv e-prints, July 

2011. 

[4] Steven R. White. Density matrix formulation for quantum renormalization groups. Phys. 

Rev. Lett, 69(19):2863-2866, Nov 1992. 
[5] U. Schollwock. The density-matrix renormalization group. Rev. Mod. Phys., 77(1):259-315, 

Apr 2005. 

[6] M. J. Hartmann, J. Prior, S. R. Clark, and M. B. Plenio. Density Matrix Renormalization 
Group in the Heisenberg Picture. Physical Review Letters, 102(5) :057202, February 2009. 
arXiv:0808.0666. 

[7] Stellan Ostlund and Stefan Rommer. Thermodynamic limit of density matrix renormalization. 

Phys. Rev. Lett, 75(19):3537-3540, Nov 1995. 
[8] G. Vidal. Efficient classical simulation of slightly entangled quantum computations. Phys. 

Rev. Lett, 91:147902, 2003. 
[9] G. Vidal. Efficient simulation of one-dimensional quantum many-body systems. Phys. Rev. 

Lett., 93:040502, 2004. 

[10] T. Barthel, C. Pineda, and J. Eisert. Contraction of fermionic operator circuits and the 
simulation of strongly correlated fermions. Phys. Rev. A, 80(4):042333, October 2009. 
arXiv:0907.3689. 

[11] R. N. C. Pfeifer, P. Corboz, O. Buerschaper, M. Aguado, M. Troyer, and G. Vidal. Simulation 
of anyons with tensor network algorithms. Phys. Rev. B, 82(11):115126, September 2010. 
arXiv: 1006.3532. 

[12] R. Konig and E. Bilgin. Anyonic entanglement renormalization. Phys. Rev. B, 82(12):125118, 
September 2010. arXiv:1006.2478. 



36 



[13] Jacob Biamonte. Lectures on Tensor Network States, QIC 890/891 Selected Advanced Top- 
ics in Quantum Information. IQC, The University of Waterloo, Waterloo Ontario, Canada 

(youtube lecture series), 2011. |http : //www . qubit . o rg/iqc2011, 
[14] A. J. Daley, S. R. Clark, D. Jaksch, and P. Zoller. Numerical analysis of coherent many-body 

currents in a single atom transistor. Phys. Rev. A, 72(4):043618, October 2005. 
[15] R. Steinigeweg, S. Langer, F. Heidrich-Meisner, I. P. McCulloch, and W. Brenig. Coherent 

spin-current oscillations in transverse magnetic fields. ArXiv e-prints, October 2010. 
[16] S. R. Clark, J. Prior, M. J. Hartmann, D. Jaksch, and M. B. Plenio. Exact matrix product 

solutions in the Heisenberg picture of an open quantum spin chain. New Journal of Physics, 

12(2):025005, February 2010. arXiv:0907.5582. 
[17] Martin Bruderer, Alexander Klein, Stephen R. Clark, and Dieter Jaksch. Polaron physics in 

optical lattices. Phys. Rev. A, 76(1):011605, Jul 2007. 
[18] M. Bruderer, T. H. Johnson, S. R. Clark, D. Jaksch, A. Posazhennikova, and W. Belzig. 

Phonon resonances in atomic currents through bose-fermi mixtures in optical lattices. Phys. 

Rev. A, 82(4):043617, Oct 2010. 
[19] J Schachenmayer, G Pupillo, and A J Daley. Time-dependent currents of one-dimensional 

bosons in an optical lattice. New Journal of Physics, 12(2):025014, 2010. 
[20] T. H. Johnson, S. R. Clark, and D. Jaksch. Dynamical simulations of classical stochastic 

systems using matrix product states. Phys. Rev. E, 82(3):036702, Sep 2010. 
[21] F. Verstraete, V. Murg, and J. I. Cirac. Matrix product states, projected entangled pair 

states, and variational renormalization group methods for quantum spin systems. Advances 

in Physics, 57:143-224, 2008. 
[22] C. V. Kraus, N. Schuch, F. Verstraete, and J. I. Cirac. Fermionic projected entangled pair 

states. Phys. Rev. A, 81(5):052338, 2010. arXiv:0904.4667. 
[23] P. Corboz, R. Orus, B. Bauer, and G. Vidal. Simulation of strongly correlated fermions 

in two spatial dimensions with fermionic projected entangled-pair states. Phys. Rev. B, 

81(16):165104, 2010. arXiv:0912.0646. 
[24] G. Vidal. Class of Quantum Many-Body States That Can Be Efficiently Simulated. Physical 

Review Letters, 101(11):110501, 2008. 
[25] G. Vidal. Entanglement Renormalization: an introduction. ArXiv e-prints, 2009. 0912.1651, 
[26] P. Corboz and G. Vidal. Fermionic multiscale entanglement renormalization ansatz. Phys. 

Rev. B, 80(16)465129, October 2009. arXiv:0907.3184. 
[27] P. Corboz, G. Evenbly, F. Verstraete, and G. Vidal. Simulation of interacting fermions with 

entanglement renormalization. Phys. Rev. A, 81(1):010303, January 2010. arXiv:0904.4151. 
[28] J. C. Baez and A. Lauda. A Prehistory of n-Categorical Physics. ArXiv e-prints, 2009. 

larXiv: 0908. 24691 

[29] A. Joyal and R. Street. The geometry of tensor calculus i. Advances in Mathematics, 88(55), 
1991. 

[30] P. Selinger. A survey of graphical languages for monoidal categories. ArXiv e-prints, August 
2009. arXiv:0908.3347. 

[31] Yves Lafont. Penrose diagrams and 2-dimensional rewriting. In Applications of Categories 
in Computer Science, volume 177 of London Mathematical Society Lecture Note Series, pages 
191-201. Cambridge University Press, 1992. 

[32] Yves Lafont. Equational reasoning with 2-dimensional diagrams. In Term Rewriting, volume 
909 of Lecture Notes in Computer Science, pages 170-195. Springer- Verlag, 1995. 

[33] Roger Penrose. Applications of negative dimensional tensors. Combinatorial Mathematics and 



37 



its Applications, Academic Press, 1971. 
[34] J. C. Baez and M. Stay. Physics, Topology, Logic and Computation Rosetta Stone. ArXiv 

e-prints, March 2009. larXiv: 0903. 03401 
[35] Yves Lafont. Towards an algebraic theory of boolean circuits. Journal of Pure and Applied 

Algebra, 184:2003, 2003. 

[36] Y. Guiraud. The three dimensions of proofs. ArXiv Mathematics e-prints, December 2006. 
|arXiv:math/0612089[ 

[37] Carolyn Brown and Graham Hutton. Categories, Allegories, and Circuit Design. In Proceed- 
ings of the 10th Annual IEEE Symposium on Logic in Computer Science. IEEE Computer 
Society Press, Los Alamitos, California, July 1994. 

[38] M. Aulbach, D. Markham, and M. Murao. The maximally entangled symmetric state in terms 
of the geometric measure. New Journal of Physics, 12(7):073025, July 2010. arXiv:1003.5643. 

[39] J. I. Cirac and F. Verstraete. Renormalization and tensor product states in spin chains and 
lattices. J. Phys. A Math. Gen., 42:4004, 2009. 

[40] J. Smith and M. Mosca. Algorithms for Quantum Computers. ArXiv e-prints, January 2010. 
arXiv: 1001. 0767. 

[41] I. Wegener. The complexity of boolean functions. Wiley- Teubner, 1987. online 

at: |http : //eccc . hpi-web . de/, 
[42] E. Boros and P.L. Hammer. Pseudo-boolean optimization. Discrete Applied Mathematics, 

123(1-3): 155-225, 2002. 

[43] Grzegorz Malinowski. Many-valued logics. Clarendon Press: Oxford University Press, 1993. 
SERIES : Oxford logic guides. 

[44] A. Kitaev, A. Shen, and M. Vyalyi. Classical and quantum computation. AMS, Graduate 
Studies in Mathematics, 47, 2002. 

[45] Michael Nielsen and Isaac Chuang. Quantum computation and quantum information. Cam- 
bridge University Press, 2000. 

[46] L. Amico, R. Fazio, A. Osterloh, and V. Vedral. Entanglement in many-body systems. Reviews 
of Modern Physics, 80:517-576, April 2008. 

[47] R. B. Griffiths, S. Wu, L. Yu, and S. M. Cohen. Atemporal diagrams for quantum circuits. 
Phys. Rev. A, 73(5):052309, May 2006. 

[48] D. Gross, J. Eisert, N. Schuch, and D. Perez-Garcia. Measurement-based quantum computa- 
tion beyond the one-way model. Phys. Rev. A, 76(5):052315, November 2007. 

[49] L. Lamata, J. Leon, D. Perez-Garcia, D. Salgado, and E. Solano. Sequential Implementa- 
tion of Global Quantum Operations. Physical Review Letters, 101(18):180506, October 2008. 
arXiv:0711.3652. 

[50] C. Schon, K. Hammerer, M. M. Wolf, J. I. Cirac, and E. Solano. Sequential generation of 
matrix-product states in cavity QED. Phys. Rev. A, 75(3):032311, March 2007. arXiv:quant- 
ph/0612101. 

[51] S. Ostlund and S. Rommer. Thermodynamic limit of density matrix renormalization. Phys. 

Rev. Lett, 75:3537, 1995. 
[52] M. Fannes, B. Nachtergaele, and R. F. Werner. Finitely correlated states on quantum spin 

chains. Lett. Math. Phys., 25:249, 1992. 
[53] F. Verstraete, V. Murg, and J. I. Cirac. Matrix product states, projected entangled pair 

states, and variational renormalization group methods for quantum spin systems. Advances 

in Physics, 57(2): 143-224, 2008. 
[54] G. Vidal. Entanglement renormalization. Phys. Rev. Lett., 99:220405, 2007. 



38 



[55] Y.-Y. Shi, L.-M. Duan, and G. Vidal. Classical simulation of quantum many-body systems 

with a tree tensor network. Phys. Rev. A, 74(2):022320, Aug 2006. 
[56] L. Tagliacozzo, G. Evenbly, and G. Vidal. Simulation of two-dimensional quantum systems 

using a tree tensor network that exploits the entropic area law. Phys. Rev. B, 80(23) :235127, 

December 2009. 

[57] Norbert Schuch, Michael M. Wolf, Prank Verstraete, and J. Ignacio Cirac. Simulation of 
quantum many-body systems with strings of operators and monte carlo tensor contractions. 
Phys. Rev. Lett, 100(4) :040501, Jan 2008. 

[58] F. Mezzacapo, N. Schuch, M. Boninsegni, and J. I. Cirac. Ground-state properties of quantum 
many-body systems: entangled-plaquette states and variational Monte Carlo. New Journal of 
Physics, 11(8):083026, August 2009. 

[59] H. J. Changlani, J. M. Kinder, C. J. Umrigar, and G. K.-L. Chan. Approximating strongly cor- 
related wave functions with correlator product states. Phys. Rev. B, 80(24) :245116, December 
2009. 

[60] John Baez et al. Bimonoids from biproducts. The n-Category Cafe Blog. online 
at: http : / /golem . ph . utexas . edu/ category/2010/ 09/bimonoids Jrom_biproducts . html. 

[61] Samson Abramsky and Bob Coecke. Categorical quantum mechanics. Chapter in the Handbook 
of Quantum Logic and Quantum Structures vol II, Elsevier, 2008. 

[62] Peter Selinger. Dagger compact closed categories and completely positive maps: (extended 
abstract). Electronic Notes in Theoretical Computer Science, 170:139 - 163, 2007. Proceedings 
of the 3rd International Workshop on Quantum Programming Languages (QPL 2005). 

[63] Discrete and switching functions. McGraw-Hill Int. Book Co., 1978. 

[64] M. Cohn. Inconsistent canonical forms of switching functions. IRE Transactions of Electronic 
Computers, 1962. 

[65] A. Mukhopadhyay and G. Schmitz. Minimization of exclusive-or and logical-equivalence 
switching circuits. IEEE Trans, on Computers, 1970. 

[66] V. Bergholm and J. D. Biamonte. Categorical quantum circuits. Journal of Physics A Math- 
ematical General, 44(24):245304, June 2011. 

[67] D. Aharonov. A simple proof that toffoli and hadamard are quantum universal. 2003. 
lquant-ph/0301040 

[68] Joachim Kock. Frobenius algebras and 2-d topological quantum field theories. Cambridge 
University Press, 2003. 

[69] P. W. Shor. Fault-tolerant quantum computation. ArXiv Quantum Physics e-prints, 1996. 

arXiv: quant-ph /960501 1 . 
[70] Eric Dennis. Toward fault-tolerant quantum computation without concatenation. Phys. Rev. 

A, 63(5):052314, Apr 2001. 
[71] A. Carboni and R.F.C. Walters. Cartesian bicategories i. Journal of Pure and Applied Algebra, 

49:11-32, 1987. 

[72] Bob Coecke and Ross Duncan. Interacting quantum observables: Categorical algebra and 

diagrammatics. aXriv preprint 0906.4725, 2011. 
[73] C. Kassel. Quantum groups. Springer Graduate Texts in Mathematics, 1994. 
[74] Samson Abramsky and Bob Coecke. A categorical semantics of quantum protocols. Proceedings 

of the 19th IEEE conference on Logic in Computer Science (LiCS'04), 2004. 
[75] J. D. Biamonte. Nonperturbative k -body to two-body commuting conversion Hamiltonians 

and embedding problem instances into Ising spins. Phys. Rev. A, 77(5):052331, May 2008. 
[76] Dhiraj K. Pradhan. A Theory of Galois Switching Functions. IEEE Tr. Computers, 27(3) :239- 



39 



248, 1978. 

[77] S.L. Hurst. Multiple- Valued Logic - Its Status and its Future. IEEE Tr. Computers, 
33(12):1160-1179, 1984. 



